General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



+ • 



(NASA-TM-73281) ANALYSIS CF GEOS-3 N79-236 15 

ALTIHETEF DATA AND EXTBACT ION CF OCEAN HAVE 

HEIG8T AND DOMINANT fc’AVELENGTE Final Bcport 

(NASA) 30 p HC A03/MF A0 1 CSCL 08C Onclas 


G 3/4 8 25202 



V. 

ik 



NASA Technical Memorandum 73281 


Analysis of GEOS-3 Altimeter Data and 
Extraction of Ocean Wave Height and Dominant 
Wavelength 


Edward J. Walsh 
Wallops Flight Center 
Wallops Island, Virginia 


NASA 

National Aeronautics and 
Space Administration 

Wallops Flight Center 

Wallops Island, Virginia 23337 


ANALYSIS OF GEOS-3 ALTIMETER DATA AND EXTRACTION OF 
OCEAN WAVE HEIGHT AND DOMINANT WAVELENGTH 


Edward J. Walsh 
NASA Wallops Flight Center 


ABSTRACT 


The GEOS-3 satellite radar altimeter has the capability to provide accurate, low 
noise measurements of the significant wave height (SWH) of ocean waves through their 
effect on the shape of the return pulse. When the amplitude and timing biases are removed 
from the GEOS-3 Sample and Hold (S&H) gates, the mean return waveforms can be excellently 
fitted with a theoretical template which represents the convolution of the radar point 
target response; the range noise (jitter) in the altimeter tracking loop; the sea surface 
height distribution; and the antenna pattern as a function of the range to mean sea level. 
Several techniques of varying complexity to remove the effect of the tracking loop jitter 
in computing the wave height are considered. They include realigning the S&H gates to 
their actual positions with respect to mean sea level before averaging; using the observed 
standard deviation on the altitude measurement to remove the integrated effect of the 
tracking loop jitter, and using a look-up table to correct for the expected value of range 
noise. Analysis of skewness in the GEOS return waveform demonstrates the potential of a 
satellite radar altimeter to determine the dominant wavelength of ocean waves. When the 
jitter in the GEOS range tracker is large, a significant reduction in range noise reduc- 
tion can be achieved through reprocessing the data. But when the tracker noise is small 
only marginal improvement is possible. 


INTRODUCTION 


The Geodynamics Experimental Ocean Satellite (GEOS-3) was launched on March 30, 1975 
into an orbit whose mean altitude is 843 km. The satellite contains a radar altimeter 
operating at 13.9 GHz which transmits 100 pulses per second towards the subsatellite 
point. The transmitted pulse of 12.5 ns nominal duration illuminates a pulse-limited spot 
on the earth's surface whose diameter is 3.6 km on flat terrain. The return waveform is 
square-law detected and point-sampled by 16 Sample and Hold (S&H) gates of 12.5 ns width 


and spaced 6.25 ns apart. One of the major goals of the satellite was to demonstrate that 
ocean wave heights along the GEOS-3 groundtrack could be extracted from the radar altimeter 
data. The potential of pulse-limited radar altimeters to measure ocean wave height has 
been discussed In detail (Barrick, 1972; Miller and Hayne, 1972; Berger, 1972; Walsh, 

1974; Brooks and Dooley, 1975), and will not be repeated here. By comparison with in-sltu 
data from waverlder buoys and a shlpborne wave recorder Walsh et a1_. (1978) have demon- 
strated that an airborne pulse-limited altimeter at 3 km altitude can accurately measure 
wave height. 

The author and a number of other GEOS sea state investigators have collaborated on a 
paper (Fedor et al_. , 1979) which Intercompares our various algorithms, determining the 
standard deviations and biases of the algorithms using both simulated and actual GEOS 
data. That paper also compares the GEOS-derlved wave heights with data buoy records. 

Those Intercomparisons show that all the algorithms are in good agreement and that the 
wave height values being produced from GEOS data are accurate. Parsons (1979) has demon- 
strated on a global basis that GEOS wave height measurements are superior to several other 
data sources utilized operationally In ship routing and weather forecasting. The purposes 
of this paper are to examine some GEOS data in detail to demonstrate that the contributions 
to the mean return waveform are well understood, to examine the techniques for extracting 
the wave height and skewness of the sea surface height distribution from it, to demonstrate 
that the skewness may be used to determine the dominant wavelength of the ocean waves, and 
to gain an appreciation for the promise of the improvements contained In the SEASAT-1 
altimeter launched on June 26, 1978. 

RADAR ALTIMETER RETURN PULSE SHAPE 


The radar mean return pulse Is the convolution of the radar system point target 
response, the jitter (range noise) in the altimeter tracking loop, the sea surface height 
distribution, and the antenna pattern as a function of the range to mean sea Isvel. 
Significant wave height (SWH) 1$ defined as the average of the heights of the one-third 
highest waves in a long sequence of waves observed at a point. For the large spatial 
averages considered in this paper (70 km) the sea surface height distribution for even a 
very narrow spectrum would still be nearly Gaussian In which case SWH can be assumed to 
equal four times the standard deviation of the height distribution (Neumann and Pierson, 
1966). To make it easier to assess their relative Importance, the widths of the various 
contributions to the radar mean return signal will be dealt with in terms of effective 
significant wave heights. For example, four times the standard deviation in meters of the 
jitter In the range tracking loop will be called the jitter significant wave height, SWH^, 
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Four times the standard deviation of the rise time of the leading edge of the GEOS-3 
return pulse (multiplied by half the speed of light to convert It to radar range) will be 
called the raw significant wave height, SWH R . If the antenna pattern Is deconvolved (for 
GEOS-3 this Is approximately equal to differentiating the leading edge of the pulse) and 
If the other contributors to the return pulse are normally distributed then 

SWH r 2 = SWHp 2 + SWHj 2 + SWH 2 (1) 

where the subscript P Identifies the effective significant wave height associated with the 
radar point target response (transmitted pulse width and any receiver broadening). 


S&H GATE AMPLITUDE AND TIMING BIASES 


The 16 GEOS S&H gates have amplitude and timing biases (Walsh, 1979) which should be 
corrected for to obtain the most accurate wave height results. Figure 1 shows histograms 
of the outputs of five (1, 6, 8, 10, 16) of the 16 S&H gates. The abscissa quantization 



Figure 1. Histograms of the outputs In mv (proportional to received power) 

of five (1, 6, 8, 10, 16) of the 16 GEOS Sample-and-Hold (S&H) gates. 
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of 7.5 mv Is approximately equal to the quantization of the analog-to-dlgltal (A/D) con- 
verters to which the signal voltages (proportional to power) of the S&H gates were applied. 
Since the ordinate Is logarithmic, an exponential distribution of power would result In a 
straight line. That Is the case for the gates sampling the noise ahead of the return 
pulse (gate 1). 

In the region where the return signal Is emerging from the noise (gates 6 and 8 In 
Figure 1) the distribution initially decays exponentially but then decreases more slowly 
for higher values. This breaking away of the tall probably accounts for the abrupt 
Increase of the fourth order moment as the signal comes out of the noise which has been 
reported by C, L. Rufenach (1976 Fall AGU Meeting) and suggested as a sensitive indication 
of the presence of a return signal. When the signal level Is high the peak of the distri- 
bution is away from the origin but the decay Is exponential. 

The A/D converter quantization varied somewhat from gate to gate. If the quantiza- 
tion was larger than 7.5 mv then, occasionally, a bln in the histogram might be empty or 
one quantization level might be split between two adjacent bins. This behavior Is In 
evidence In the histogram for gate 16 in the bins around 120, 210, 262.5 and 315 mv. 

Histograms of the outputs of the first five S&H gates are shown in Figure 2 using a 
much smaller (0.1 mv) bin size. In general, the observations appear as spikes spaced at 
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Figure 2. Histograms of the outputs of the first five S&H gates using 
a much smaller bin size (0.1 mv) than Figure 1. 


roughly 7.5 mv Intervals. Some of the spikes are narrow (gate 4) while some are rela- 
tlvely broad and exhibit structure of their own (gate 1). The substructure would not have 
an Important effect on the data since even the broad spikes are narrow with respect to the 
7.5 mv spacing. The Important consideration Is that the various gates have significant 
biases with respect to each other. 

The lowest entry In the histogram for gate 1 Is 5 mv lower than the lowest entry for 
gate 2. If the populations of the corresponding spikes In the histograms of gates 1 and 2 
are approximately the same then that shift In the distribution would result In the same 
magnitude difference In the mean values. Looking at the value of the first spike In each 
histogram was the first cut at determining amplitude biases In the gates. In the early 
gates a refinement could be made If the populations of the spikes were not the same, such 
as gates 4 and 5. If an exponential process were being sampled one could Interpolate to 
find where the process ended. For gate 5 that would have been 4.8 mv below gate 4 rather 
than the 6 mv that would be Indicated by the lowest spike In gate 5. 

The estimates of the biases were refined using an Iterative technique to be described 
later. The results are shown In the top of Figure 3 where the average of the biases for 
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Figure 3. Amplitude and timing biases for the GEOS altimeter S&H gates. 
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all gates has been arbitrarily set to zero. The bottom of Figure 3 shows the timing 
biases for the various gates, also zero meaned, which were determined by E. L. Hofmelster 
of the General Electric Company from prelaunch test data (private coranunicatlon, 1976). 
Hofmelster had reported earlier (GEOS- 3 PI Meeting/1975) that gate 13 was reading low and 
was highly correlated with gate 12 and uncorrelated with gate 14. Both effects are 
explained by the timing bias. The top of Figure 4 shows data points (dots) for three 
successive 960-pulse averages of GEOS data and the theoretical returns (solid lines) 
resulting from matching the data with a template formed by convolving a skewed normal 
distribution (representing the combined effects of surface height distribution, radar 
point target response, and range noise) with an exponential decay representing the antenna 
pattern. The amplitude biases are most apparent In the saw tooth pattern of the first 
five gates. 


PASS 4594 

5.99 M 5.54 M 6.24 M 



4.78 M 5.22 M 5.40 M 



Figure 4. Data points (dots) for three successive 960-pulse averages of 
GEOS data and theoretical returns (solid lines) resulting from 
curve fitting before and after the biases of Figure 3 were 
removed. 

The bottom of Figure 4 shows the same data with the biases removed and the fit redone. 
The most important timing bias is that gate 13 Is closer to gate 12 than It is to gate 14. 
The fit of the data to the theoretical curves is much improved and the estimates of SWH R 
have all lowered. Figure 5 is similar to Figure 4 except that the sea state was higher. 
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PASS 4523 

11.88 1* 10.35 M 9.87 M 



10.60 M 9.82 M 9.16 M 



Figure 5. Data points (dots) for three successive 960-pulse averages of 
GEOS data and theoretical returns (solid lines) resulting from 
curve fitting before and after the biases of Figure 3 were 
removed. 


The GEOS-3 data is recorded In frames containing 320 pulses. Three-frame averages 
were used in the data analysis so that the standard deviation In the radar data due to 
Rayleigh fluctuations would be only 3.2 percent of the mean value at each point. Refine- 
ments to the amplitude biases were made by applying the timing biases and the preliminary 
amplitude biases to approximately 30 three-frame averages and determining the template. 

If the amplitude biases were correct the data points should fluctuate about the theoreti- 
cal returns in a random manner with zero mean. If the mean difference between the data 
points and the curves is not zero the amplitude biases were adjusted accordingly and the 
procedure repeated. It is apparent from the fit of the data in Figures 4 and 5 that 
removing the biases In the GEOS data should produce more accurate results. The timing and 
amplitude biases were removed from all the data considered in this paper. 


PROCEDURES FOR FITTING DATA WITH THEORETICAL TEMPLATE 


The parameters of the theoretical template for the return pulse were varied In an 
iterative manner and compared with the data to minimize a fit parameter (Walsh et al_. , 
1978). If the fit parameter Is the sum of the squares of the residuals of the data to the 
template then a least-squares (LS) fit Is being performed. If the residuals are nor- 
malized by the standard deviation observed at each data point before squaring then the fit 
Is a maximum likelihood estimate (MLE). Both procedures were used In fitting the GEOS 
data. Little difference was found between the two methods for high sea states. At low 
wave heights the MLE sometimes had trouble converging because relatively small residual 
errors In the amplitude biases In the noise gates would be highly weighted due to the 
small standard deviations In those gates. Then the solution would stop before It con- 
verged to the correct wave height because the fit parameter would Indicate that not much 
Improvement was being made. The biases and the curves In Figures 4 and 5 were derived 
using MLE estimates but LS estimates were employed In the rest of the paper. J. F. R. 
Gower (1979) has Indicated that the variance of the SWH estimates can be reduced using a 
hybrid fit where the residuals are normalized by the standard deviation but a limit Is 
placed on how small the standard deviation may become. 


TECHNIQUES FOR COMPENSATING FOR RANGE NOISE 


Figure 6 shows schematically the Ideal and actual positioning of the S&H gates by the 
tracker. Ideally the gates should be slewed In range smoothly so that gate 10 would 
always be positioned at the range corresponding to mean sea level (MSL). Actually the 
tracker does not adjust the range smoothly but jitters the gates back and forth about 
their Ideal positions (shown exaggerated In the figure) so that they are somewhat mis- 
aligned. This range noise is generally assumed to be normally distributed. In practice, 
however, the deviation of the GEOS S&H gates from their Ideal position is generally ckewed 
and sometimes can even differ significantly from a skewed-normal distribution. 

In the high data rate mode of GEOS when all of the Individual pulse-return values of 
the S&H gates and the range servo error (rse) are available the range jitter effect can be 
eliminated entirely. Each rse tells the actual change in range of the S&H gates from 
their previous position and the rso can be integrated to determine the exact spatial 
positioning of the S&H gates for each pjlse. Over a short distance the GEOS-3 satellite 
altitude variation is very nearly linear. If the Integrated rse is detrended by per- 
forming a LS fit to t straight line then the residuals are the mispositioning of the gates 
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Figure 6. Schematic representation of Ideal and actual positioning 
of GEOS S&H gates by the range tracker. 

Indicated on the right side of Figure 6. This procedure also removes any linear trend in 
MSL due to surface slopes but It should be emphasized that the detrending does not destroy 
any information. Any sea surface slope that Is present can be recovered If the actual 
satellite range rate Is known. 

The detrending was done over groups of 960 pulses which constitutes three frames of 
GEOS data and covers approximately a 70 km ground track. Using only one frame {320 pulses) 
was Insufficient to Insure that the anomalously large deviations (to be discussed later) 
sometimes observed In the tracker would not adversely affect the results. The residuals 
of the detrended integral of the rse can be used to eliminate range tracker noise by 
realigning the individual pulse returns before they are averaged. One can envision on the 
right side of Figure 6 that the rse residuals allow the S&H gates to be relabeled to a 
primed set of gates that as a fixed with respect to MSL. The data in the pr1„ied reference 
frame would then be free of range tracker noise. If, for example, the gates were posi- 
tioned so that S&H gate 10 was actually at the position gate 5 should have been with 
respect to MSL, then the data coming In thrcugh gate 10 for that pulse would be placed In 
the bin corresponding to gate 5 prime for the tt/erage data. That essentially turns the 
actual tracker shown at the right In Figure 6 Into the ideal tracker shown at the left and 
eliminates the SWHj term in (1). For realigned data (1) reduces to 

SWH R 2 « SWH p 2 + SWH 2 (realigned data). (2) 
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Equation (2) Is the equation of a circle, as is (1) If the effects of pulse and 
jitter are combined in a single term such as 

SWripj 2 - $WH p Z f SWHj 2 . (3) 

o 

If both sides of (2) are divided by SWH R the circle has unit radius as shown In Figure 7. 
It Is apparent that when the SWH Is large compared to the pulse {or pulse and jitter in 
the case of (1)) It is essentially equal to SWH R . In that region the curve Is nearly 
horizontal so that the error In estimating SWH Is essentially equal to the error In 
estimating SWH R . 



SWHp/SWH R 

Figure 7. Graphical representation of the relationship between SWH, 

SWHp, and SWH r . 

Even In perfectly realigned data the effects of Rayleigh fading In the radar return 
pulse will cause the measured SWH R to differ from the Ideal. The standard deviation of 
that measurement depends on a number of things and It Is inversely proportional to the 
square root of the number of pulses averaged. Even a small percentage error In estimating 
$WH R produces a large percentage error In SWH when SWH Is small compared to the pulse 
since the curve In Figure 7 Is nearly vertical. Figure 7 shows that only a five percent 
over estimate In SWH R when the SWH was actually zero would cause the calculated SWH to 
equal 0.31 SWH p (1.2 m for GEOS). Similarly, If SWH R were under estimated by five percent 
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for zero SWH then (2) Indicates that the SWH Is Imaginary and equal In magnitude to 0. 31 p. 
It Is apparent that the altimeter pulse width Is fundamental In determining the noise In 
the measurement of low values of SWH. Since the SEASAT-1 altimeter pulse width Is 
nominally a factor of four smaller than the GEOS pulse. Its noise level when measuring low 
wave heights will be correspondingly Improved. 

Having the SWH computation result In an Imaginary number should produce no consterna- 
tion and Introduce no error when one realizes that It Is simply a mathematical artifact. 
Over a series of measurements the measured value of SWH R will tend to be high as often as 
It Is low, so making the Imaginary values of SWH negative Is one way of dealing with the 
problem. Then a series of values of SWH could be filtered to obtain « better estimate of 
the actual SWH and the positive and negative values of SWH wou^ tend to cancel If the 
actual SWH were zero. But since the transformations from SWH R to SWH using (1) or (2) are 
highly non-linear when SWH Is near zero the best procedure to use In filtering the data 
would be to filter the values of SWH R to obtain the bast estimate of that quantity and 
tien make the transformation to SWH, setting It equal to zero If the computed value were 
inary. 


DETERMINATION OF SWH 


Figure 8 shows dashed curves representing the raw significant wave height, SWH R , and 
the range jitter, SWH JP obtained from non-overlapping three-frame averages of GEOS data. 
The effective significant wave height of the radar point target response, SWH p , is con- 
stant at 3.81 m and will not be shown on the figures. Computing SWH using (1) results In 
the data points (squares and diamonds) which Indicate a fairly constant SWH around 3.5 m 

except for the average of frames 67, 68, 69 where It is indicated as zero. At th-,t point 
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(1) indicated that SWH was equal to -13.9. In this Instance the imaginary SWH .alue is 
not due to Rayleigh fluctuations in the return pulse but is caused by over estimating SWHj 
because the range jitter Is not normally distributed In that region. It can be eliminated 
by realigning the individual return pulses before averaging so SWH Is determined from (2) 
instead of trying to compensate for the integrated effect of the tracker misalignment in 
the averaged data using (1). The SWH values determined using the realigned data are shown 
by the solid curve in Figure 8. The circle is the SWH measured by a 3 ns pulse-limited 
radar altimeter (Hughes Aircraft Company, 1976) from a NASA WFC C-54 aircraft during a 
simultaneous GEOS underflight at 3 km altitude. The two methods of extracting SWH from 
GEOS data are in reasonable agreement except where the SWH value was anomalously low using 
the first technique. 

The range noise for the 12 frames of data corresponding to the four diamonds in 
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Figure 8 Is shown at the top of Figure 9. The residuals of the detrended Integral of the 
rse are shown at the top of the figure. At the bottom of the figure are shown histograms 
of the deviations of the S&H gate positions from the ideal position for each of the four 
three-frame averages. The abscissa and ordinate of each histogram have been normalized 
appropriately using the standard deviation observed in the data to correspond to the 
normal distributions indicated by the dots. The first and second histograms are very 
nearly normally distributed with some skewness and the fourth histogram is not too far 
off. But the third distribution Is non-normal and considerably more narrow than is 
indicated by its standard deviation. The reason for the anomalously low SWH value in 
Figure 8 was because the high value of SWHj was not a true indication of the width of the 
distribution. The large excursion is an example of the tracker anomalies referred to 
earlier and is not caused by a variation in the level of the sea surface. 

Because a relatively small number of large residuals can cause too large a value of 
SWHj there is a tendency for the values of SWH produced by (1) to be lower than those 
obtained using (2) even when the effect is not striking as it was in Figure 8. Figure 10 
shows a pass where no imaginery SWH values were produced by (1) but that algorithm tended 
to produce lower values than the realignment process. The higher values produced by (2) 
are in better agreement with the value of SWH measured from 3 km altitude by the Naval 
Research Laboratory radar altimeter (also on the NASA aircraft, see Walsh et aK (1978) 
for an analysis of its measurement capabilities). Figure 11 shows the range noise and 
histograms for the four three-frame averages indicated by the diamonds on Figure 10. The 
second of the four histograms is the closest to a skewed-normal distribution and yet the 
results of (1) and (2) differ by 1.5 m. 

Figure 12 shows a third case where there is both good and poor agreement between the 
two techniques. Since SWH is a measure of how much energy is contained in the sea surface 
wave pattern it generally would be expected to be a smoothly varying function when aver- 
aged over the many ocean wavelengths that would be contained in the 70 km distance traveled 
during a three-frame average. One would expect that most of the fluctuation in SWH R from 
one three-frame average to another would be caused by, and in phase with, fluctuations in 
range noise. This behavior is in evidence in all three of the figures, 8, 10, and 12, but 
a particularly good example appears in the Figure 12 in the region from frame ^2 to frame 
45. There is a saw-tooth variation in both SWH R and SWHj but the SWH determined from (1) 
shows a gradual trend. It is not surprising that the agreement between SWH determined by 
(1) and (2) is good in that region. 

On the other hand, there is also a saw-tooth pattern in evidence in the region from 
frame 61 to frame 75. Although the SWH R and SWHj patterns are in phase, the high values 
of SWHj are anomalously high due to those distributions of rse residuals being non-normal 
as evidenced by the histograms in Figure 13. The result is that the SWH calculated from 
(1) decreased sharply when SWH R increased at frames 64 and 70. 
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SWH (M) 


PASS 4381 



Figure 8. Plots of SWH R , SWHj and SWH for non-overlapping three-frame 
averages of GEOS data for pass 4381. 
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Figure 9. A plot of the residuals of the detrended integral of the 
range servo error and histograms of the residuals for the 
data corresponding to the diamonds in Figure 8. 




Figure 10. Plots of SWH^, SWHj and SWH for non-overlapping three-frame 
averages of GEOS data for pass 4523. 



Figure 11. A plot of the residuals of the detrended integral of the 
range servo error and histograms of the residuals for the 
data corresponding to the diamonds in Figure 10. 
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Figure 12. Plots Of SWHp, SWHj and SWH for non-overlapping three-frame 
averages of GEOS data for pass 4481. 



Figure 13. A plot of the residuals of the detrended integral of the 
range servo error and histograms of the residuals for the 
data corresponding to the diamonds in Figure 12. 


It Is also unrealistic from a consideration of the hardware that SWH would decrease 
sharply when SWH^ Increases. The range tracker bandwidth increases as the SWH decreases 
and the faster response should decrease the range noise. A decrease in the actual SWH 
should decrease both SWHj and SWH R , not Increase them. The range noise distributions at 
frames 67 and 73 are closer to being normal and the agreement Is better between the values 
of SWH determined from (1) and (2). 

During the preprocessing of the GEOS data by NASA WFC the altitude standard devia- 
tion, a^, is computed for each frame and placed on the data tapes before they are dis- 
seminated. Figure 14 is a scatter plot, versus SWHj, of the altitude standard deviation 
multiplied by four to convert it to an effective wave height. Because SWHj was calculated 
on a three-frame basis, the square root of the mean of the squares of the three values 
was used for the comparison. The slope of the straight line through the data points is 
45° but the SWHj intercept is 0.5 m which indicates that the two quantities typically 
differ by 0.5 m. If all the values of rse were not available (low data rate) or if one 
didn't wish to take the computer time to integrate the rse and detrend it, the values of 
o R could be used to determine a value of SWHj for use in (1). 



Figure 14. Scatter plot of altitude standard deviation computed by 
NASA WFC during preprocessing of GEOS data versus SWHj. 
The straight line is at a 45° angle. 
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Hofmeister et al_. (1976) have done a detailed analysis of the GEOS range tracker. 

But a very simple analysis will suffice here to put the GEOS data In perspective. To 
first order one would expect the tracking jitter to be Inversely proportional to the mean 
slope of the leading edge of the return pulse. Another way of saying the same thing is 
that it should be proportional to the rise time of the return pulse. But that Implies 


SWHj « K ^SWHp 2 + SWH 2 (4) 

where K is a constant which is equal to the ratio SWHj/SWH p when SWH is zero. 

Figure 15 is a scatter plot of range jitter versus the SWH determined from realigned 
data using (2). The open circles indicate points associated with tracker anomalies. The 
dashed curve is a plot of (4) for K = .6955. Figure 16 is a scatter plot of SWH deter- 
mined from realigned data versus SWH p from non-aligned data. The dashed curve is a plot 
the equation resulting from substituting (4) into (1). 

SWH R 2 = (1 + K 2 ) (SWHp 2 + SWH 2 ) (5) 


The curve is in fairly good agreement with the data points but a better fit to the data 
could be made using a curve whose slope was slightly steeper at the high end. 



o I I i 1 i - J l 1 1 1 1 

01 8348078 # 10 

SWH FROM MALMNCD DATA (M) 

Figure 15. Scatter plot of range jitter versus SWH determined from realigned 
data. The open circles indicate points associated with tracker 
anomalies. The dashed curve is a theoretical result assuming 
tracker noise is proportional to rise time of return pulse. 
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The most accurate technique of extracting SWH from GEOS-3 data would be to do the 
realignment and use (2). When that approach is not possible (low data rate) or when the 
computer time to do the realignment is too costly then the best procedure would probably 
be to use a look-up table to go directly from SWH R to SWH rather than use (1). The look- 
up table could be generated using a curve similar to the one shown in Figure 16 but one 
which had been least-squares fitted to the data points. Then SWHj would not have to be 
considered. Even regions of non-normally distributed range noise would cause little 
difficulty since they affect SWHj much more than they do SWH^. 



SWH r TROM NONALIGNED DATA (M) 

Figure 16. Scatter plot of SWH determined from realigned data versus 

SWHp from nonaligned data. The dashed curve is the theoretical 
result expected using the dashed curve from Figure 15. 


18 


Figures 8, 10, and 12 show values of SWHj reaching 5 to 7 m. Under similar condi- 
tions the value of SWHj for the SEASAT altimeter should not exceed 0.40 m. That, combined 
with the low value of SWH p for SEASAT, indicates that SWH will essentially equal SWH R for 
all but the lowest wave heights. Ai;d they should readily be determined from a look-up 
table. 


MEASUREMENT OF SKEWNESS OF SEA SURFACE HEIGHT DISTRIBUTION 


As discussed earlier, the theoretical template which was fit to the GEOS data con- 
tained skewness as a parameter. The skewness determined from the template would be the 
raw skewness, X R , associated with SWH R . Prof. W. J. Pierson has shown (private communi 
cation, 1977) that If each of the contributors to SWH R is sxev/ed- normal distribution 
then A R is equal to 



SWHj 


( 6 ) 


where X, A p , and Aj are the skewness values associated with the surface height distribu- 
tion, the radar point target response, and the range noise. 

As is apparent in the histogram of Figure 9, 11, and 13, the GEOS range noise is 
generally highly skewed and sometimes non-normally distributed. The mean value of Aj is 
0.5 and the standard deviation is 0.4. Figure 15 shows that SWHj is roughly the same 
magnitude as SWH for high wave heights and considerably larger than SWH for low wave 
heights. So range noise would make it very difficult to use (6) to extract A from GEOS 
data. If realigned data is used then (6) reduces to (7). 


= A 


SWH 



(7) 


Figure 17 shows some results of a simulation performed by the author and described in 
Fedor e£ aK (1978). Rayleigh fluctuations were added to theoretical GEOS mean return 
wave forms for realigned data of known SWH R and X R simulate returns corresponding to 
averages of various numbers of pulses. Estimates of the various parameters were made 
using the simulated data and their means and standard deviations were computed. The 
standard deviation of X R was approximately independent of the mean value and varied 
inversely as the square root of the number of pulses averaged. Figure 17 shows that 
estimation of skewness is a noisy process for one-frame averages of GEOS data, especially 
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Figure 17. Expected standard deviation of A R estimate derived from 
stimulated realigned GEOS data. 


for the lower values of SWH. 

To determine the on-orbit value of skewness for the radar point target response 
GEOS-3 data was processed for three passes (4590, 5443, 6629) which same passed by closest 
approach {< 30 km) a NOAA National Data Buoy (EB-15) Indicating that the wave field was 
low swell (wind < 2.4 m/s, SWH < 0.7 m, wave period > 5 seconds). The GEOS data itself 
indicated that there was a large area in the vicinity of the buoy on each of the passes 
where the wave height was very low. 

The mean value of A R for the three-frame averages of realigned GEOS data was 0.46 for 

the segments used from the three passes. The standard deviation on those X R estimates was 

0.26 which agrees with the simulation results since it is just V3 times smaller than the 

value indicated In Figure 17 for the 0 to 1 m range of SWH for one-frame averages. 

The theory discussed in the next section indicates that A would be approximately 0.1 
or less for the conditions described by the buoy. The small ratio of SWH/SWH R and the 
cubic weighting in (7) make the contribution of A to A R negligible under those conditions 
and the value of A p can be calculated to be 0.48. 

Equation (7) can be rearranged to form 

/SWH R \ 3 / SWHp 

x = x r I swfT ) ‘ x p i swR“ 



When SWH is large compared tc SWH p then it Is approximately equal to A R . But when SWH is 
small compared to SWH p then A is, in general, the difference of two relatively large terms 
and would be difficult to determine accurately. 
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DETERMINATION OF DOMINANT WAVELENGTH 


In the recent theoretical study quantitatively confirmed by laboratory and field 
data, Huang and Long (1979) determined that the skewness of the height distribution of 
waves Is proportional to the significant slope, 5. This significant slope was defined as 
the ratio of rms wave height to dominant wavelength where the dominant wavelength was 
determined by the maximum of the wave height frequency spectrum. Their theoretical result 
was 


X * 8ir §, (9) 

Data from the passes shown In Figures 10 and 12 provide an excellent opportunity to 
examine the possibility of using their result to determine dominant wavelength. The wave 
height was essentially all sea In Figure 10 at the position of the ground truth data point 
and all swell at the same position In Figure 12. 

For the position of the NASA ground truth aircraft along pass 4523 (Figure 10) the 
hindcast Indicated that the wind was 20 m/s at 240° and the wind-driven sea had an 8 m SWH 
with an eight second period while the swell was only 1 m with a 14 second period. The 
hindcast SWH seems to be accurate since the composite SWH would be 8 m which Is in agree- 
ment with both the GEOS SWH and the aircraft measured SWH. Two Independent estimates of 
the wind were obtained from the aircraft (Robert Mennella, private communication, \S1). 
Taking 0.8 times the wind velocity measured by the aircraft at 150 m altitude resulted In 
a surface wind estimate of 22 m/s at 229°, Estimating the magnitude of the surface wind 
using the NRL radar in the Wind-Wave-Radar (WWR) mode (Hammond et al_. , 1977) resulted in 
18 m/s. 

For the position of the aircraft along pass 4481 (Figure 12) the hindcast indicated 
that the wind was 5 m/s at 260° and the wind-driven sea was only 0.5 m with three second 
period while the swell was 6 m with ten second period. The hindcast also seems fairly 
accurate In this case since the GEOS and aircraft measured SWH values were 5.4 m. The 
wind estimate using 0.8 times the aircraft-measured wind velocity at 150 m resulted 1r: 

9 m/s at 281° and the WWR estimate of the surface wind magnitude was 7.6 m/s. Since Che 
highest SWH that could be expected to be produced by an 8 m/s wind would still be less 
than 1.5 m it is reasonable to assume that the waves were essentially all swell. 

Figure 18 is a plot of X for passes 4523 and 4481 versus distance from the ground 
truth data points. The values of X were determined from nine-frame averages of GEOS data 
(210 km ground track) to Improve the accuracy of the measurement. The center of the nine- 
frame average was shifted three frames at a time (70 km) so that adjacent data points on 
the plot were obtained from average data which had six frames in conmon. Both plots show 
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DISTANCE FROM GROUND TRUTH (KM) 

Figure 18. Variation of sea surface skewness as determined from GEOS 
data for wind-driven sea and swell conditions. The mean 
values of X over the data points were 0.48 for the wind- 
driven sea and 0.21 for the swell. 

a significant variation In skewness, but the mean of the values for the wind-driven sea 
(pass 4523) Is 0.54 which Is significantly higher than the 0.21 mean value for the swell 
(pass 4481). 

Figure 19 Is a plot of (7) where two axes have been used on the abscissa. One axis 
Is in terms of S and one Is In terms of the ratio of dominant wavelength to SWH, which Is 
0.25 times the reciprocal of §. Data from passes 4523 and 4481 are shown on the figure 
with the vertical positions of the diamond and circle Indicating the mean of the respec- 
tive skewness values shown In Figure 18. The error bars on the mean values are the 
standard deviations of the Figure 18 data points reduced by the square root of one third 
the number of points averaged to account for the correlation caused by the overlap In the 
averages. The horizontal position of the diamond Is the ratio of the wavelength indicated 
by the hlndcast (96 m wavelength corresponding to an eight second period) to the 8 m SWH 
determined by GEOS (which agreed with the hindcast and aircraft data). The horizontal 
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position of the circle uses the ISO m wavelength corresponding to the ten second hlndcast 
period and the 5.4 m GEOS SWH (which agreed with the aircraft data but was slightly lower 
than the 6 m hlndcast SWH). 

Figure 19 Indicates that an altimeter could be used to determine the dominant wave- 
length of the ocean waves. The SWH and skewness can be determined from the returned pulse 
shape. Then SWH can be multiplied by the ratio of dominant wavelength to SWH determined 
from Figure 19 and the skewness. Two data points are by no means definitive but they do 
demonstrate the potential. An extensive analysis of satellite radar altimeter data and 
comparison with surface truth needs to be done to firmly establish this remote sensing 
technique. The pulse repetition rate of the SEASAT-1 altimeter is ten times higher than 
GEOS-3. the pulse width Is four times smaller* the data point spacing on the return 
waveform Is twice as dense and the SWHj Is below 0.4 m. Therefore, for SEASAT-1 X should 
nearly equal X R , the spatial sampling of X on the ocean will be higher and the noise level 
on the measurement will be lower. 


NOISE REDUCTION IN RANGE MEASUREMENT 


In addition to SWH R and X R , another parameter that Is determined In the process of 
fitting the theoretical return to the data Is the arrival time (epoch) of the return 
pulse. The left side of Figure 20 shows that when the GEOS raw range noise Is high then 
realigning the pulses and determining the epoch for 30-pulse averages will significantly 
reduce tf* noise In the range measurement compared to simply averaging the range data. 
However, if the GEOS tracker noise Is In Its typical range, then the amount of Improvement 
to be gained using the more complex technique Is small (right side of Figure 20). This 
same conclusion was arrived at by Dooley et ah (1979) after a more complex analysis. It 
should be noted that range tracker anomalies such as the one seen in Figure 9 completely 
disappear when the pulses are realigned and retracked. That is the best evidence that 
they are not indications of surface elevation changes. Realigning pulses should never 
be necessary on SEASAT-1 altimeter data since the range noise should generally be below 
10 cm. 


CONCLUSIONS 


The contributions to the mean return waveform of a pulse-limited radar altimeter are 
well understood. The GEOS-3 satellite radar altimeter Is providing accurate, low noise 
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Figure 20. GEOS tracker range noise and range noise resulting from 

averaging range data points versus realigning and retracking 
the return pulses. 


measurements of the significant wave height of ocean waves. Analysis of skewness In the 
GEOS return waveform demonstrates the potential of a satellite radar altimeter to deter- 
mine the dominant wavelength of ocean waves. When the jitter in the GEOS range tracker Is 
large, a significant reduction In range noise reduction can be achieved through reproces- 
sing the data. But when the tracker noise Is small only marginal Improvement Is possible. 
The SEASAT-1 altimeter should produce even more accurate SWH and dominant wavelength 
values with greater spatial resolution than GE0S-3. 
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